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Abstract 

We present a detailed analysis of the hard X-ray (>2.5 
keV) EPIC spectra from the first observations of NGC 
4151 made by XMM-Newton. We fit the spectra with 
a model consisting of a power-law continuum modified 
by linc -of-sight absorption (arising in both partially pho- 



in a substantial (Nh ~ 10 



23 cm 2 ) line-of-sight gas col- 



umn density (e.g. 


Holt ct al. 1980 




Yaqoob et al. 1993 


Weaver et al. 1994a 


3 


Weaver et al. 


19941). 



toioniz bd and neutral gas) plus additional iron-K emission 
and absorption features. This model provides an excellent 
overall fit to the EPIC spectra. The iron Ka line is well 
modelled as a narrow Gaussian component. In contrast 
to several earlier studies based on data from ASCA, a 
relativistically broadened iron Ka emission feature is not 
required by the XMM-Newton data. The upper limit on 
the flux contained in any additional broad line is ^8% of 
that in the narrow line. The measured intrinsic line width 
(a = 32 ± 7 eV) may be ascribed to (i) the doublet nature 
of the iron Ka line and (ii) emission from low ionization 
states of iron, ranging from neutral up to FeXVII. The 
additional iron absorption edge arises in cool material and 
implies a factor ~ 2 overabundance of iron in this compo- 
nent. 
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1. Introduction 

The Seyfert 1 galaxy NGC 4151 is one of the brightest 
Active Galactic Nuclei (hereafter AGN) accessible in the 
X-ray band, and has been extensively studied by all major 
X-ray missions. This observational focus has revealed the 
X-ray spectrum of NGC 4151 to be a complex mixture of 
emission and absorption components, originating at a va- 
riety of locations from the innermost parts of the putative 
accretion disk, out to the extended narrow-line region of 
the galaxy. The X-ray spectrum emanating from the active 
nucleus is dominated by an intrinsic X-ray to 7-ray contin- 
uum that appears to be produced by the thermal Comp- 



Some earlier studies based on ASCA observations have 
concluded that the line profile is complex, and is com- 
posed of with an intrinsically narrow component (Weave: 
1994b| ; Zdziarski et al. 1994, etc) plus a relativis 
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tically broadened line feature (e.g. Yaqoob et al. 1993 



tonization of soft seed photons (e.g. 


Haardt & Maraschi 


1991 


Zdziarski et al. 1994 


, Zdziarski ct al. 1996, 


Zdziarski 


et al. 2000; Pctrucci et al. 2000]) , plus additional contribu- 



tions from reprocessing in the form of Compton-reflection 
and features associated with neutral iron (e.g. iron-K fluo- 
rescence and a neutral iron edge) . Below ~ 5 keV the hard 
continuum is strongly cut-off by photoelectric absorption 



Wang et al. 2001). Furthermore the broad line profile 
was reported to be variable on timescales of 10 4 s, corre- 
sponding to an emitting region of <0.02 AU, suggestive of 
an origin close to the supermassive black hole presumably 
in the inner regions of a putative accretion disk. 

In this paper we focus on the iron-K features present in 
the very high signal-to-noise EPIC spectra of NGC 4151. 
Our analysis is based on the spectral "template" model, 
developed by [Bchurch fc Warwick (2002) to interpret ear- 
lier ASCA and BeppoS AX observations. 



2. The XMM-Newton observations 

NGC 4151 was observed with XMM-Newton three times 
during the 21 s * and 22 Tld of December, 2000 (orbit 190). 
NGC 4151 was p ositioned on-axis in both the EPIC MPS 
and PN cameras ( |Turner et al. 200l| ; [Btrudcr ct al. 200l|) 



with the medium filter in place. Of the three EPIC ob- 
servations, one (~33 ks) was performed with the MOS 
and PN in small window mode, whilst the remaining two 
observations (~63 and ~23 ks) were operated in full win- 
dow mode for both instruments. The recorded events were 
screened with the latest available release of the XMM- 
Newton Science Analysis Software (SAS) to remove known 
hot pixels and other data flagged as bad. The data were 
processed using the latest CCD gain values and only X-ray 
events corresponding to patterns 0-12 in the MOS cam- 
eras and 0-4 in the PN camera were accepted. Investiga- 
tion of the full field count-rate for all three observations 
revealed a single flaring event during the first observation. 
This flaring event was screened from the data, resulting in 
effective total exposure times of 110 ks and 91 ks respec- 
tively for the MOS and PN instruments. An investigation 
into the impact of pile-up on the observation showed that 
the effect was negligible in all the observations, largely 
due to the relatively faint nature of NGC 4151 during this 
period. 
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A source and background lightcurve were extracted 
from each individual, screened observation. These lightcurves 
revealed that NGC 4151 remained remarkably constant 
(to within 10%) for the duration of the observations, al- 
lowing us to analyse the XMM-Newton data without be- 
ing overly concerned with the considerable spectral vari- 
ability that characterized the previous ASCA "long- look" 
observation ( ^church fc Warwick 2002 ). Source and back- 
ground spectra were extracted from the same regions as 
the lightcurves. The spectra from the two full window 
mode observations were co- added to produce a single source 
spectrum and a single background spectrum for each in- 
strument. The spectra were binned to a minimum of 20 
counts per spectral channel, in order to apply x 2 minimi- 
sation techniques. Here we consider only the 2.5-12 keV 
spectrum. 



Table 1. Best- fit spectral model parameters 
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3. The spectral analysis 
3.1. The initial spectral template model 



1 



We adopt the spect ral template model described by, Schurch 
fc Warwick (2002) which includes the following emission 
components: 

A power-law continuum with a normalization, A\ , and 
photon index, T, exhibiting a high-energy break at 100 
keV; 

A neutral Compton-refl ection component (modelled by 
PEXRAV in XSPEC, [Magdziarz fc Zdziarski 1995|) 



with only the reflection scaling factor, R, as a free pa- 
rameter. The parameters relating to the incident con- 
tinuum were tied to those of the hard power-law com- 
ponent. In addition cos i was fixed at 0.5 and the metal 
abundance in the reflector was fixed at the solar value; 
3. An iron Ka emission line of intensity Jx a at an energy 
Eku with an intrinsic line width ctkq- 

The complex absorber is represented as the product of 
two absorption components, namely a warm column den- 
sity A^H,warm and a cold gas column iVn,coid- For further 



details of the photoionization modelling see Schurch & 
Warwick (20"02j| and priffiths et al. (1998)| . The adopted 



spectral model includes absorption arising in the line-of- 
sight column density through our own Galaxy, applied to 
all three emission components (iVH,Gai= 2xl0 20 cm -2 ). 



3.2. The spectral fitting 

The template spectral model was fit simultaneously to 
both the small window mode and the co-added, large win- 
dow mode EPIC spectra. The initial model allowed the 
continuum and reflection normalizations, the line proper- 
ties (centroid energy, width and normalization) and the 
ionization parameter of the warm absorbing column to be 
free parameters. The values for all the remaining model 
parameters were set at those used in Bchurch fc War 



wick (2002) . The spectra, along with the spectral template 



Energy (keV) 

Figure 1. The XMM-Newton EPIC spectra of NGC 4151 and the 
spectral template model. The MOS data, shown here in black (MOS 
1) and red (MOS 2), are from the small window mode observation. 
The PN data, shown here in green, is from the Full Window mode 
co-added observation. 

model fit, is shown Figure [l] and the best-fit parameters 
are given m Table 0, Model 1. The quoted errors (here as 
elsewhere in this paper) are at the 90% confidence level 
as defined by a Ax 2 =2.71 criterion (i.e., assuming one in- 
teresting parameter). This model gave an acceptable fit 
to the spectrum (% 2 =5465 for 4973 d.o.f), although the 
data/model residuals suggest an over- prediction of the flux 
between 7 and 8 keV. A closer inspection of the shape of 
the data in this region reveals a sharp drop at ^7 keV, 
which we interpret as an absorption edge from neutral 
iron. Including an additional edge in our model, with the 
edge energy and the optical depth a free parameters, im- 
proved the fit dramatically (Ax 2 = 173 for 2 d.o.f; Figure 
^|) . The best-fit parameters for the model incorporating an 
additional edge are given in Table |, Model 2. 

Including a further contribution to the iron Ka line 
profile from a broad (<r=0.2 keV) component does not 
improve the model fit, and sets an upper limit on the flux 
in the broad component of l.OxlO -5 photons cm~ 2 s" 1 . 
Replacing the Gaussian line in the best-fit model with the 
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Figure 2. The XMM-Newton EPIC spectra of NGC 4151 fit with 
the modified spectral template model. The EPIC spectra, MOS 1 
(black), — MOS 2 (red) and PN (green), are fit with the Hftectral tern. 



ment of the fit and results in an upper limit on the flux 
in the broad component of ~8% of that in the narrow 
component. Modelling the line profile with a relativisti- 
cally broadened disk-line profile with parameters values 
as quoted in |Yaqoob et al. (1995) results in a poor fit 
to the data (x =8800 for 4971 d.o.f). Similarly, modelling 
the line profile with the parameters quoted by Wang et 
al. (2001) also results in a poor fit to the data (x 2 =9438 
for 4971 d.o.f). In fact the broad line profiles proposed 
by Wang et al. and Yaqoob et al. are both ruled out at 
>99.99% significance. 

The measured flux in the narrow iron Ka line is 1.26^o 03 > 
photons cm~ 2 s _1 , significantly lower (at ^> 4tr) than 
that measured in the ASCA, BeppoSAX and Ginga ob- 



servations where the line flux was more typically 2.2_ 



plate <r>u,c QjXdjfiomjjX prigi feature £mm neaxtxal agoa n * Z 1 

keV. 



+0.2 , 
OA. 



10 4 photons cm 2 s 1 ; (Schurch fc Warwick 2002, Yaqoob] 



Warwick 1991) the present measurment is also some- 



what lower than the line flux reported in more recent 
Chandra observations (1.8xl0 -4 photons cm -2 s _1 ; Ogle 



broad emission from an accretion disk (disk-line model 
in XSPEC, Fabian et al. 1989] ) results in a similar fit to 
that for a narrow line (x 2 =5309 for 4970 d.o.f). However, 
in this case the derived inner radius of the accretion disk, 
r;, ~1000r g gives a disk-line profile that closely resembles 
the observed narrow line. 



4. Discussion 

The hard X-ray spectrum of NGC 4151 is well modelled by 
the spectral template model developed in Schurch fc War 



wick (2002) , with surprisingly few modifications. In the 
XMM-Newton observations NGC 4151 was in a relatively 
faint state as judged against previous extensive monitoring 
of the source by missions such as EXOSA T and Ginga. The 
continuum is ~60% weaker, than the faintest level mea- 
sured during earlier BeppoSAX and ASCA observations, 
whereas the level of the Compton reflection component is 
~60% stronger. The ionization parameter of the warm ab- 
sorber falls within the range of ionization states recorded 
during the ASCA long-look observation (£=2.48-2.65). 

4.f. The nature of the iron line 



et al. 2000), implying a possible decrease in the line flux 
over a period of ^1 year. 

XMM-Newton measures a line width of 32± 7 7 cV, some- 
what greater than the line width measured by Chandra (^5 
eV) . Possible explanations for the non-zero intrinsic width 
of the Gaussian profile include: (i) the signature of the 
doublet nature of the underlying iron Ka line and (ii) 
emission from low ionization states of iron, ranging up 
to FeXVII, present in the warm absorber. However, mod- 
elling the line in XMM-Newton with two intrinsically nar- 
row components to represent the iron Ka doublet yields a 
marginally worse fit to the data (x 2 =5320 for 4970 d.o.f). 
The lines were fixed at the correct rest energies of the Kai 
(6.392 keV) and Ka 2 (6.405 keV) lines and the branch- 
ing ratio was fixed at the theoretical value of 1:2 respec- 
tively. We note that Ogle et al. (2000)| find no evidence 
for Doppler broadening from their spectra. 

4.2. The iron edge - a clue to iron abundance. 



Following the analysis described in Schurch & Warwick 



(2002), we model the iron Ka line in the XMM-Newton spec- 



tra with a simple Gaussian model at an energy correspond- 
ing to fluorescence of cold, neutral iron. The data are well 
modelled by this line profile and we find no evidence for 
the presence of a relativistically broadened iron Ka line. 



This is in contrast with some previous claims (Yaqoob et 
al. 1995; Wang et al. 2001) based on the analysis of 



ASCA observations. However, we note that further anal- 



ysis of these earlier data by Schurch & Warwick (2002) 



do not require the presence of a broad line. Including a 
broad line in our modelling gives no statistical improve- 



The spectral model developed in Bchurch fc Warwick (2002) 
includes iron edge features from both the neutral and 
warm absorbers and from the Compton reflection of neu- 
tral gas, assuming in each case solar metal abundances. 
These solar abundance edges are insufficient to model the 
7-8 keV region of the XMM-Newton spectrum, prompt- 
ing the inclusion of an additional edge. The derived edge 
energy is consistant with absorption by neutral iron. This 
feature can be interpreted either as a result of abundance 
gradients or as a result of a general overabundance of 
iron in the absorbing material. Previous analyses of NGC 
4151, albeit using a partial covering model, have also re- 
ported evidence for an overabundance of iron in NGC 
4151. Yaqoob et al. (I993)| quote a canonical value of 
Ap e ~2.5 times solar, based on the average of many mea- 
surements. Interpreting the additional edge required by 
the XMM-Newton data as an overabundance of iron in 
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the neutral absorber (Morrison & McCammon 1983), the 
optical depth of the edge leads to a value for the iron abun- 
dance, Ape, of 3.6j^'g 5 times solar. Alternativly, allowing 
the iron abundance of the Compton reflection component 
to replace the additional edge in the model results in an 
more modest iron overabundance (Af c ~2.3), consistent 
with the results of Yaqoob ct al. (1993)| , however in this 
case the fit becomes somewhat worse above 8 keV. If we al- 
low an overabundance in both the neutral absorption and 
in the neutral reflection component, a value of Ap c ~2 is 
required. 

5. The Conclusions 

Remarkably, the hard X-ray (2.5-10 keV) EPIC spectra 
from the recent XMM-Newton observation of NGC 4151 
are extremely well modelled by the model developed in 



Schurch fc Warwick (2002) with only minor modifications. 



The success of this simple model in modelling both previ- 
ous and current data supports that the spectral variabil- 
ity exhibited by NGC 4151 is the product of rather subtle 
changes, specifically in the ionisation state of the majority 
of the line of sight warm absorber. The parameters of the 
warm absorber fit to the XMM-Newton data are consis- 
tent with the results from the spectral analysis of the past 
seven years of ASCA and BeppoSAX observations. 

The chief modification to the template spectral model 



( Schurch & Warwick 2002 ) is the addition of an edge fea- 
ture from neutral iron. Interpreting the edge as the result 
of an overabundance of iron in the cold absorber and cold 
reflector, then a value of Ap ^2 is required. 

This analysis finds no requirement for a relativistically 
broadened iron Ka line feature in the XMM-Newton data. 
The upper limit on the possible flux present in any addi- 
tional broad line is ~8% of that in the narrow line compo- 
nent. There is evidence that the line flux in NGC 4151 has 
decreased recently, possibly due to the fact that the source 
has been fainter during the last few years than was pre- 
viously observed. The low luminosity (Lx,2-io ~2xl0 42 
erg s _1 ) of the current observation is consistent with this 
view. 

The intrinsic non-zero width of the best-fit Gaussian 
line profile is suggestive of a further level of complexity to 
the iron Ka line feature, with possible explanations for the 
non-zero intrinsic width including: (i) the signature of the 
doublet nature of the iron Ka line, (ii) emission from low 
ionization states of iron, ranging up to FeXVII, present 
in the warm absorber. 
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